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Abstract. 

We present a simple technique for studying collisions of ultracold atoms in the 
presence of a magnetic field and radio- frequency radiation (rf). Resonant control of 
scattering properties can be achieved by using rf to couple a colliding pair of atoms 
to a bound state. We show, using the example of ^Li, that in some ranges of rf 
frequency and magnetic field this can be done without giving rise to losses. We also 
show that halo molecules of large spatial extent require much less rf power than deeply 
bound states. Another way to exert resonant control is with a set of rf-coupled bound 
states, linked to the colliding pair through the molecular interactions that give rise 
to magnetically tunable Feshbach resonances. This was recently demonstrated for 
^^Rb [T]. We examine the underlying atomic and molecular physics which made this 
possible. Lastly, we consider the control that may be exerted over atomic collisions by 
placing atoms in superpositions of Zeeman states, and suggest that it could be useful 
where small changes in scattering length are required. We suggest other species for 
which rf and magnetic field control could together provide a useful tuning mechanism. 
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1. Introduction 

Many recent studies of ultracold gases liave depended on tlie manipulation of atomic 
collisions with external fields. One way of achieving such control is magnetically tunable 
Feshbach resonances [21 |3l H] . Also, lasers have been used to create optical Feshbach 
resonances by coupling a colliding pair to a bound state of an excited potential [5l Ej . 
A number of recent works have used optical frequency lasers in combination with 
magnetically tunable Feshbach resonances to probe [71 [8] and modify [9] scattering 
and bound state properties. This can provide additional control capabilities, due 
to advantageous properties of laser beams such as the ready tunability of power 
and frequency. Furthermore, some experiments of current interest, such as color 
superfluidity [10] and Efimov physics [11], could benefit from the ability to independently 
control scattering lengths between different pairs of a multicomponent gas. 

Radiofrequency radiation (rf) is an essential and much used tool in atomic physics, 
being central to the operation of experiments such as atomic clocks [12] . In the context 
of ultracold gases, rf has been used both to dissociate [13] and associate [T^ [T5l [16] 
molecules, as well as to drive transitions between bound states [17|. As a probe, rf allows 
measurements of interaction effects [I8l[l9], molecular binding energies [131 [HI |15] , and 
the pairing gap of fermionic superfluids [20]. These experiments creating or studying 
bound states used the rf to create a degeneracy between dressed bound and scattering 
states, i.e., a Feshbach resonance. This provides a clear motivation for study of how the 
same degeneracy can be used to control the collision properties of a pair of atoms. 

In Ref . [1] , rf control of ®^Rb collisions was demonstrated and studied theoretically. 
Furthermore, rf has been considered as a means for controlling scattering lengths in a 
variety of scenarios, with a number of different theoretical methods. Moerdijk et al. 
considered collisions of rf-dressed Na atoms in a magnetic trap using a coupled-channels 
approach [21]. They showed that favourable conditions for evaporative cooling could 
be created, and suggested controlling scattering lengths using rf coupling to create a 
Feshbach resonance. Zhang et al. used a two-channel parametrisation of a magnetically 
tunable Feshbach resonance, together with atomic rf dressing to tune the resonance 
location, or multi-frequency rf to independently tune different scattering lengths in a 
multicomponent gas [22]- Tscherbul et al. performed a coupled channels analysis of 
^^Rb collisions, studying both the creation and manipulation of resonances with rf [23] . 
Papoular et al. suggested rf control of collisions in gases at zero magnetic field |24j . 
Lastly, we note the suggestions of Alyabyshev et al. to use rf to control atom-molecule 
interactions [25], particularly to suppress inelastic collisions [26] . 

Studies of ultracold gases have benefited immensely from theoretical prediction and 
characterisation of Feshbach resonances [27] , for which coupled channels calculations are 
a leading technique [2H1 121 [22]. These have the drawback of being computationally 
intense in some cases. For this reason, a number of simpler techniques have been 
developed, including the asymptotic bound state [30l [31] , accumulated phase [32] , and 
three-parameter van der Waals models [33]. Our technique [33] is able to calculate 
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scattering and bound state properties on the basis of only three parameters describing 
the interactions. In Ref. [1] we extended it to include the effects of rf radiation on 
collisions of ®^Rb. In this paper we present several new results illustrating the utility of 
rf for controlling ultracold collisions, and explain our model in greater detail. 

Radiofrequency radiation can modify collision properties in three ways, each of 
which we examine here. Firstly, rf can couple a colliding pair and a bound state 
(bound-free coupling), as used in the above-mentioned experiments on creation and 
measurement of molecules fi3\ [111 IISI [IE]. This creates a situation similar to an 
optical Feshbach resonance. We consider nonzero magnetic fields, and large halo states 
created by magnetically tunable Feshbach resonances, giving example results for ^Li. 
We show that this reduces the power required for useful control from that necessary for 
a deeply bound state. In general, rf resonances cause losses by coupling the colliding 
pair to a large number of energetically lower configurations of the internal atomic states, 
which we refer to as channels. Interestingly, though, we have found some regions of rf 
frequency and magnetic field in which creating a resonance using a deeply bound state 
does not lead to losses. Secondly, rf can couple several bound states together (bound- 
bound coupling), which can then interact with the colliding pair through the molecular 
interactions which give rise to magnetically tunable Feshbach resonances. This was 
recently demonstrated for ^^Rb in Ref. [1]. The Franck-Condon factors between two 
bound states are generally larger than those between bound and scattering states, 
allowing rf control with comparatively low power. Our theory was applied to this 
situation in Ref [1]; here, we give a more detailed analysis. Lastly, dressing an atom with 
rf creates a superposition of atomic states. The collision of two such atoms then involves 
a superposition of several entrance channels (free-free coupling) . If the scattering length 
varies significantly between channels, this alone could be used to vary the scattering 
length. 

To make this paper self-contained, we start with a description of our three- 
parameter model in Sec. [2l The main extension necessary to include the effects of rf is 
the formulation of an appropriate basis for the scattering calculation, which we discuss 
in Sec. [31 Our results for bound-free, bound-bound and free-free coupling are presented 
in Sees. 14.11 14.21 and 14.31 respectively. For each, we consider the extent of control that 
is possible, the rf power required, and the accompanying losses. We conclude in Sec. [5l 

2. Scattering theory approach 

We consider the collision of two ^S alkali atoms. In the presence of a magnetic field, 
B = Bz, the atoms are described by a set of Zeeman states |a) with energies Ea- We 
label such states a = a,b,. . ., in order of increasing energy. These correlate with the 
zero field states \f,mf), where / is the total angular momentum of the atom, and ruf 
its projection along z. At non-zero field, mj remains a good quantum number but / 
does not. The collision of a pair of atoms is described by an expansion of the interaction 
Hamiltonian in channels |ai + Q;2) \(^iTii), defined by the state of atoms 1 and 2, the partial 
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wave i of their collision, and m^, the projection of i along z. Several coupled channels 
are typically involved in a collision. The threshold energy of a channel is defined as 
Eai+a2 = + Ea^, i-G. the energy of two atoms in the relevant atomic states at 
asymptotically large separation r, with zero relative kinetic energy. If the threshold 
energy of a channel is below the total energy of the colliding atoms, the channel is 
called open; otherwise, it is referred to as closed. At zero magnetic field, channels of the 
same total angular momentum T = fi + f2 + i are coupled. At non-zero field, however, 
channels of the same Mt = mi + m2 + are coupled. 

The valence electrons of the two colliding atoms give rise to two Born-Oppenheimer 
(BO) interaction potentials, one each of singlet and triplet symmetry. These are 
isotropic, independent of £ and m^, and for large r take the form of a van der Waals 
potential, —Cq/t^. Here, Cq is the van der Waals coefficient, which is the same for both 
BO potentials. Expanding the Hamiltonian in terms of the channels \ai + 02) the 
diagonal matrix elements have the long range form 

where irir is the reduced mass. Off-diagonal elements of the BO potentials decay 
exponentially as r increases. These can lead to inelastic spin relaxation (ISR), the loss 
of atoms by decay into an energetically lower channel. Such losses are also referred to 
as inelastic spin-exchange collisions. Weaker effects such as relativistic spin- dependent 
interactions are not included in the present calculations. The potential takes the 
form of Eq. ([1]) for r > r*, where r* is a distance at which the splittings between diagonal 
elements have the same order of magnitude as the off-diagonal elements. This distance 
is usually of order 20 ao, where = 0.05292 nm is the Bohr radius. 

For r < r*, the depths of and splittings between the singlet and triplet BO potentials 
are much larger than the relative kinetic energy of the colliding atoms and the atomic 
hyperfine splittings. This energy-scale separation leads to the concept of a quantum 
defect, in which a simple parametrisation is used to account for the short-range physics, 
and is matched with solutions of the long-range potential. These long-range solutions 
are much easier to obtain. Such an approach is valid because the short range region 
is not probed in detail by the colliding pair. We note that quantum defect theory has 
been used in a wide variety of contexts including atomic scattering [Ml EH] , electron-ion 
collisions and Rydberg states [361 [37], and nucleon scattering [38] . 

Our approach is based on a series of papers by Gao (see Refs. [391 HDl |35l E] and 
references therein), in which a number of powerful tools for pure C„/r" potentials were 
presented. For each channel we define a reference potential, taken to be Eq. (IT|) extended 
over all r. From this we calculate reference functions / and g, which are two linearly 
independent solutions of the Schrodinger equation. We use these to form vectors / and 
g spanning all channels. In the original problem, mixing between channels occurred due 
to spin exchange. Here, this mixing is incorporated in the r = boundary condition. 
As done in the work of Gao [35], by choosing / and g appropriately we can express this 
boundary condition as a short range K matrix, K'-*\ that is independent of collision 
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energy E and partial wave. Here, we use bold font to indicate a matrix. In fact, the 
multichannel wavefunction with our approximate Hamiltonian can be written as 

^{r) = f{r) - K(^)^(r) , (2) 

for all r. 

A convenient molecular basis for calculating K*^'*^ is that described by 
\{siS2)S{iii2)I] Fi; TMt), where Si_2 and ii^2 are the electronic and nuclear spin angular 
momentum of atoms 1 and 2, respectively. The two electron spins are coupled together, 
as are the two nuclear spins, to give the total electron spin S = S1 + S2 and total nuclear 
spin I = ii + 12- These are then coupled to give F = S + I , which is coupled to the 
partial wave £ to give the total angular momentum, T. In this basis, K^'^^ is diagonal, 
with diagonal entries depending only on whether the corresponding channel is of singlet 
(S* = 0) or triplet (5 = 1) symmetry. Their values i^s,t are given by the relation [35] 

/- /TrA^s.t + tan(7r/8) 

«s,t/a = V2— — — — . (3) 

As,t - tan(7r/8) 

Here, a^^t are the scattering lengths of the singlet and triplet potentials, respectively, 
a = 2-i/2[r(3/4)/r(5/4)]/^dw is the mean scattering length, Z^dw = {2m,CQ/hY^'^/2 
is the van der Waals length, and T{z) is the gamma function. We use a frame 
transformation [l2l |33] to convert the K matrix to the basis \ai + a2)\i'm£). As these 
are not associated with singlet or triplet symmetry, we then have off-diagonal terms in 

Determination of scattering properties is based on the calculation of the reference 
functions / and g. The Schrodinger equation with a pure van der Waals potential can 
be solved analytically [35]. At short range, / and g have the limiting form 

f(s) = -^s'/'cos{V2s-'-n/A), (4) 

g{s) = - -^s'/' sin(v^s-2 - 7r/4) . (5) 

Here, we have defined s = r//vdw- These functions can be recognised as the zero-energy 
Wigner- Kramer- Brillouin solutions to the van der Waals potential. For open channels, 
they are linked to the long range forms. 



f{r) ^ 



9{r] 



Zfs{E) sin ( A;r - j - Zfc{E) cos ikr- — 
Zgs{E) sin ( /cr - y j - Zgc{E) cos ( A;r - y 



(6) 



(7) 



Here, k is the wavenumber of the relative motion. We note that the Z coefficients 
are also i dependent. For closed channels, the reference functions become analogous 
superpositions of exponentially growing and decaying functions. 

The reference functions allow us to find the physical K(£') matrix for the open 
channels, from which scattering properties can be extracted. This takes the form [35] 

K{E) = -[Zf,{E) - Z,,(E)Keff][Z;,(E) - Zg,{E)K,s]-' , (8) 
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where 

Keff = kW + IxiE) - KW]"1kW . (9) 

Here, the 'oo' and 'cc' refer to the open and closed channel blocks of K.^'^\ while 'oc' and 
'co' indicate the open-closed blocks. The Z matrices are diagonal and have entries equal 
to the Z coefficients of Eqs. (Q and (jT]) evaluated for the corresponding channels. The 
diagonal matrix of the bound state phase in the closed channels, x{E), is analogously 
found from the long range behaviour of the reference functions in the closed channels. 
It allows the calculation of bound state energies from the determinental equation [35] 

det{x{E)-Kil^) = 0. (10) 

We are able to extract all desired scattering properties from K(£'). We first 
calculate the S matrix, S{E) = [1 + zK(£')][l — zK(£')]~^. In the presence of just 
one open channel and the limit of A; — j- 0, the scattering length can then be found from 
the relation 

S{E) = exp{-2ika) . (11) 

We note that S(-E'), which is only defined for open channels, is a scalar for this case. 
Equation ( iTTl) can still be used when there are several open channels. However, the 
diagonal S matrix element of the entrance channel will have less than unit modulus. 
We can then reinterpret the right side of Eq. ( ITT!) as exp{—2ika) , where a = a — ib is the 
complex scattering length, and a and b are real. For any collision energy, the two-body 
decay rate coefficient is 

'<^-^,E\s.m\ (12) 

where the index i ranges over all open channels other than the entrance channel e. For 
the multichannel case, an alternative way of finding b is to extract it from K2 in the 
limit /c — )■ 0, using the formula 

Our model can be optimised to measured data by using the singlet and triplet 
scattering lengths as fit parameters [33]. This can be used for prediction of further 
resonances, or just to offset the limitations of our simplified approach. These limitations 
were discussed in Ref . [33] . Briefly, resonances arising as a result of deeply bound states 
can be predicted only approximately, as non-van der Waals parts of the potential are 
signiflcant. For bound states accurately reproduced by the —Ce/r^ potentials, however, 
our approach is accurate. The examples considered in the present work fall within this 
category. With the approximations of using a van der Waals potential and assuming 
energy independence at short range, our approach can predict all observable scattering 
properties from known atomic parameters and three properties of the interactions: Qg, 
at and Cg. 
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2.1. Decaying resonances 

Resonances invariably create some losses in ultracold collisions. Inelastic spin relaxation 
(ISR), spin-spin dipole coupling between partial waves, and three-body recombination 
are all enhanced near a Feshbach resonance. As shown in Eq. (fT3|) . two-body decay into 
energetically lower exit channels leads to an imaginary part of the scattering length. The 
losses in our calculations represent ISR. The energy gap between channels is typically 
large enough that atoms undergoing ISR are lost from the system. Our calculations do 
not include three-body recombination or decay into other partial waves. 

Theoretical work on decaying resonances has focused on the optical case, where a 
laser couples a colliding pair to an excited bound state which can spontaneously decay. 
This theory can be readily adapted to a magnetically tunable decaying resonance, or 
one in which both a magnetic field and rf are used. Following Bohn and Julienne 
we can write the complex scattering length in the limit — )■ as 

(1/2)0 - 

Here, Cbg is the background scattering length, representing the scattering length of the 
entrance channel in the absence of a resonance. We have expressed the decay rate of 
the bound state, 7, in magnetic field units, 73 = %//Xres5 where /ires is the difference 
in magnetic moment between the entrance channel and the bound state causing the 
resonance. The resonance length Ores is defined by Ores 7b = Obg^- The width and 
magnetic field location of the resonance are given by A and 5o, respectively. In the limit 
of 7 — )■ the above formulas reduce to the standard relation describing the scattering 
length around a non-decaying resonance: 

a{B) = [1 - ■ (16) 



3. rf-dressed basis 



Radio-frequency radiation drives transitions between atomic Zeeman states. In the 
two-body picture, rf couples Zeeman channels together. We therefore use a basis of 
rf-dressed channels, \ai + 02, N), where TV is the number of photons in the rf field. In 
the following sections we consider only s-wave collisions, and so omit the partial wave 
labels. In Fig. [1] we show a schematic of our coupling scheme. All channels with the 
same are coupled together by spin exchange, the molecular interactions which give 
rise to magnetically tunable Feshbach resonances and ISR in the absence of rf. We 
refer to a set of channels with the same {M^, N) as a spin-exchange block, indicated 
by a box in Fig. [H Arrows indicate rf coupling. It is usually possible to identify the 
entrance channel in which the gas and rf field are initially prepared, and for which we 
set Mt = Mo and = A'o- 
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-1 1 

Mt - Mo 



Figure 1. Schematic of our coupling sclicmc, showing photon number N vs. total 
angular momentum projection Mt- Boxes represent groups of rf-dressed Zeeman 
channels with the same {Mt,N), which we call spin-exchange blocks. The entrance 
channel is in the central {Mq,Nq) block. Arrows indicate rf coupling. A <Tx rf photon 
can couple in four additional spin-exchange blocks, with a second photon adding 
another 8. Circular polarised light, by contrast, constrains the system to one diagonal, 
as indicated by the labels. 



The polarisation of the rf field significantly alters the collisional and loss properties 
by determining which blocks couple. Light that has ax polarisation, i.e. linear and 
perpendicular to the magnetic field quantisation axis, drives (AMt = ±1, AA^ = ±1) 
transitions, as sketched in Fig. [H Here, the two '±' signs are independent - that is, 
both an increase and decrease in Mj- can be achieved with either an absorption or an 
emission of an rf photon. As a consequence, ax rf will always couple energetically lower 
exit channels, and so create losses. This can be seen from Fig. [1] if a channel |a -|- /3, N) 
is part of the basis, it is coupled to \a + f3,N — 2). When a transition from \a + f3,N) 
to \a + (3,N — 2) takes place, 2hjjj of energy is transferred from the rf field to the kinetic 
energy of the atom pair. This is typically enough energy for the atoms to be lost from the 
trap. Circularly polarised light, however, will constrain the sign of AMt with emission 
or absorption. This corresponds to moving on only one diagonal of Fig. [H which makes 
it possible to create a situation in which a pair of atoms in the entrance channel have 
no allowed exit channels. For a^ and ax light, the change in Mt with absorption or 
emission means that only spin-exchange blocks for which [Mt — Mq) + (A^ — A'^o) is even 
appear. 

We can not make a rotating wave approximation in the rf-dressed channel energies, 
because channels with off-resonant energies can contain a near-resonant bound state. 
To make sure no such bound states are omitted, we begin with the block containing the 
entrance channel and add all blocks that can be coupled by the chosen rf polarisation, 
up to the required number of rf transitions. The number of rf transitions that must be 
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included for convergence of a calculation depends on the two-body spectrum, as well as 
the strength of the rf radiation and its detuning from the atomic transition frequencies. 
The rf coupling between Zeeman channels is taken to be that between the constituent 
atomic states, including the appropriate two-body symmetrisation. The rf adds a term 
to the Hamiltonian of the form Hj-f = — (/Ii -|- /I2) ■ -Brf, where /Ii,2 are the magnetic 
moments of atoms 1 and 2, and i?rf is the rf magnetic field operator. The two-body 
matrix element can be easily formed from the matrix elements between rf-dressed atomic 
Zeeman states \a,N). For o"^ radiation, the atomic matrix element is 

{a, N\- fl- 5rf|a', N') = -Sn,n'±i^ {ne{a\Sx\a') + . (17) 

Here, is the amplitude of the oscillating field, the TV dependence of which will be 
neglected since N is assumed to be large, = (S+i—S^i) / a/2, and /ie,n are the magnetic 
moments of the electronic state and the nucleus, respectively. For cr^ polarised light, 
Sr^ and Ix must be replaced with S±i and I±i. The remaining matrix elements are 
calculated by decomposing the Zeeman states into the components of the electron spin 
s and nuclear spin i, and realising that 

{s'm'JSg\sms) = Ss'sSm'„{m.s+q)i~'^yV s(s + nis + q, -g,m,) .(18) 

Here, C is a Clebsch-Gordan coefficient, and g = —1,0, or 1 p5]. The matrix element 
for Ig can be calculated in the same manner. In this paper, we find it convenient to 
express the strength of the rf coupling in terms of the Rabi frequency between the two 
energetically lowest Zeeman states, with the pertinent polarisation: 

hn = 2\{a,No\- fl- B,i\b,No-l)\. (19) 

The Rabi frequencies for all other atomic transitions then follow from angular 
momentum algebra. This definition is weakly dependent on the bias field B, however 
its variation over the field ranges shown in our examples is negligible. 

In extending the method of Sec. [2] to include rf, we add to the approximations that 
the rf coupling is negligible at short range. That is, the matrix K^^^ is still independent 
of energy, and is also diagonal in A^. In order to use the established tools of scattering 
theory, a basis must be chosen that makes the Hamiltonian diagonal at asymptotically 
large r. For the present case, a basis of two-body Zeeman channels will have off-diagonal 
rf coupling matrix elements, as well as Zeeman energies on the diagonal. We therefore 
explicitly diagonalise the Hamiltonian for the case r — )■ 00, which provides a basis of 
rf coupled states. The eigenvectors can then be used to express the scattering matrix 
in this basis. An example of dressed ^^Rb channel energies for magnetic fields around 
8.6 G and rf frequency of 6 MHz are shown in Fig. O For this example, which we discuss 
at length in Sec. 14. 2[ the Zeeman effect is mainly linear, leading to several nearby 
avoided crossings between channels with different photon numbers. In Fig. [2] we have 
only shown the channels corresponding to the B = limit (/i = 1) + (/2 = 2). However, 
our calculations include all Zeeman channels of the relevant {Mx, N) blocks. 
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Figure 2. Radio- frequency dressed channel energies of two ®^Rb atoms as a function of 
magnetic field. Energies are given relative to i?o, the energy of the \a + e, Nq) entrance 
channel in the absence of coupling to other channels. Avoided crossings between many 
of these states occur near 8.6 G. The inset shows the same quantities over a wider 
range. Groups of channels are labelled with the number of photons associated with 
them. Here, the rf oscillation frequency is 6 MHz, the Rabi frequency is 38 kHz, and 
IG = 10-* T. 

4. Results 

The technique we have developed is general and can be used to examine any pair of 
alkali atoms. In this section we present three significantly different examples. In Sec. 14.11 
we consider directly creating resonances by rf coupling a colliding pair to a bound state, 
focussing on the example of ^Li. In Sec 14.21 we discuss ^^Rb, in which several bound 
states are coupled together by rf. This bound state manifold is then coupled to the 
colliding pair by molecular interactions. Lastly, we consider the control that can be 
obtained with rf dressing of atomic Zeeman states of ^Li in Sec. 14.31 

^Li: bound-free coupling 

The broad s-wave resonances in the a + b, a + c and b + c channels of ^Li have 
been used for studying effects such as the BEC-BCS crossover |l6l STJ HE] and Efimov 
physics lini [501 151] . Three- component ^Li gases (with atoms in Zeeman states a, b, and 
c) near these broad, overlapping resonances have also been considered as candidates 
for observing color superfiuidity [lO]. This and other applications could benefit from a 
second degree of control to allow tuning of the interactions between different component 
pairs. In this subsection we consider rf resonances in which the colliding pair of atoms 
is coupled to a bound state. The strongest bound-free Franck-Condon overlaps are 
provided by halo states. Substantial control of scattering properties is then possible 
with less power than is required for a deeply bound state. However, for some bound 
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Figure 3. Real part of the scattering length (top) and two-body loss rate cocfScient 
(bottom) for rf induced resonances in ^Li. On the left plot, these are shown as a 
function of magnetic field, with fixed rf frequencies indicated in MHz. On the right 
plot, these are shown as a function of rf frequency, with fixed magnetic fields given in 
G. Lines of the same colour correspond to the same resonance, with the five values of rf 
frequency / magnetic field in the two plots chosen to pass through the resonance peak. 
The Rabi frequency of the cr^ rf is 100 kHz, and the entrance channel is \a + b, Nq). 
The dashed line shows the effect of using rf of the same modulation amplitude. The 
resonance positions coincide with the energy of the bound state causing the resonance 
in the \a + c,Nq — 1) channel at 690 G. The width of the resonance increases as the 
bound state moves closer to threshold, enhancing the Franck-Condon overlap of bound 
and scattering states. 

states there exist ranges of rf frequency and magnetic field for which the entrance channel 
is the energetically lowest. This requires the use of circularly polarised rf, but allows 
control without the creation of losses. 

We show an example of an rf-induced resonance with a halo state in Fig. [31 This 
resonance is created in the scattering of ^Li atoms in the \a + b,NQ) channel, with rf 
coupling the colliding pair to the bound state that causes a resonance in the energetically 
closed \a + c,Nq — 1) channel at 690 G. Here, we have used a Rabi frequency of 100 kHz, 
corresponding to an oscillation amplitude of 0.5 G. The frequency of the cr~ rf is allowed 
to vary from 82 MHz to 83 MHz, which resonantly couples the bound state at magnetic 
fields in the range 600 G - 620 G. The two panels of Fig. |3] show that the resonance 
can be tuned with rf frequency or magnetic field. A larger rf frequency resonantly 
couples the bound state at a higher magnetic field, for which its energy is closer to the 
\a + c, Nq — 1) threshold. The resonance created is then wider as a function of both B 
and u, as well as providing a larger variation in scattering length and a higher peak 
loss rate. These effects are consequences of the bound state having a more halo-like 
character, and a larger bound-free overlap with scattering states of all open channels. 
These rf-coupled exit channels are the only available means of two-body decay, since in 
the absence of rf a -|- 6 is the lowest channel with Mt = 0. 

As can be seen from Eqs. ( !Ti|) and (fT5|) . the extrema of the scattering length in 
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the vicinity of such a decaying resonance are abg ± ctrcs- The resonance length Ores 
thus provides a convenient indication of how significantly the scattering length can be 
controlled with a given resonance. For uj/2'k = 82 MHz and Bq = 601.4 G we have 
Ores = 280 ao, whereas for u/2n = 83 MHz and Bq = 620.1 G we have a^^s = 13000 ao- 
The use of circularly polarised rf can reduce the number of available exit channels, and 
therefore reduce total losses. The dashed lines in Fig. |3] illustrate this for the present 
case, giving the scattering properties for rf of the same modulation amplitude as the 
solid a~ lines. However, the loss rates could still be too high for experimentally relevant 
densities. 

For some ranges of magnetic field strength and rf frequency circularly polarised 
rf allows the coupling in of a bound state without creating losses. One such example 
is shown in Fig. HI A deeply bound Mt = 1 state supported by higher channels is 
degenerate with the a + b collision threshold close to 543 G. In the vicinity of this 
magnetic field, it may be coupled to the \a + b, Nq) entrance channel using radiation, 
without opening any exit channels. This is illustrated by the dot-dashed lines in 
Fig. m The divergence in scattering length characteristic of non-decaying resonances 
is observed. In contrast to this, rf of the same frequency and power produces a 
decaying resonance, shown by the solid lines in Fig. |H A maximum in these losses 
is observed, with a background loss rate that increases as the rf frequency is brought 
closer to the atomic transition frequency. Only a slight difference in the width of the 
resonance is seen as the rf frequency is varied within the range shown, as the properties 
of the bound state are only weakly varying. This calculation uses VI/2tt = 1 MHz, 
corresponding to a modulation amplitude of 5 G. Much lower Rabi frequencies could 
be sufficient, provided that the control of the bias magnetic field is good enough. The 
frequency of the rf may be used to decide on the location of the resonance for magnetic 
field tunability, and the intensity used to choose the strength, analogous to an optical 
Feshbach resonance. 

As we discussed in Sec. [3l creating a resonance with rf always gives rise to 
losses. However, the ability to control the location and width of the loss feature could 
make such resonances useful as a knife for evaporative cooling. This was suggested for 
magnetically tunable resonances by Mathey et al. [52]. They found that the width of a 
resonance limits the temperature to which it can be used to cool a gas. While the width 
of a magnetically tunable resonance is set by molecular properties, the width of an rf 
resonance can be controlled by changing the rf power. Figure [5] shows the scattering 
length and loss rate coefficient as functions of magnetic field, for Rabi frequencies up 
to 1 MHz. Here, we have used an oscillation frequency of 77.5 MHz and (Tx polarisation. 
The Rabi frequency does not alter the nature of the resonantly coupled bound state, 
and therefore does not change the maximum variation in scattering length. For each of 
the resonances in Fig. [5] we calculate a resonance length of approximately 190 ao, close 
to the background scattering length of 184 oq. However, increasing the Rabi frequency 
increases the width of the resonance feature. For cooling, the width of the resonance 
could be reduced as the temperature is lowered. The lowest achievable temperature 
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Figure 4. Scattering length (top panel) and two-body loss rate coefficient (bottom 
panel) as a function of magnetic field, with |a-|-5, iVo) entrance channel, a Rabi 
frequency of 1 MHz and rf frequencies indicated in MHz. The rf resonances are created 
by coupling to a deeply bound state. While a.^ rf (solid lines) produces strong losses, 
as shown in the bottom panel, for (dot-dashed lines) the entrance channel is the 
energetically lowest state. Consequently, no losses are created by the rf. The small 
features at the right of the plot are due to a transition, and are therefore only 
present for the ax rf. 



would then be set by technical considerations such as magnetic field control. 

The loss properties of an optical or magnetically tunable Feshbach resonance are 
set by the decay rate of the bare bound state. This is typically almost constant in the 
vicinity of the resonance. In the optical case, this allows the stimulated coupling to 
be strengthened by increasing the laser power, increasing the resonance length. In our 
calculations, the losses are also created by stimulated coupling - i.e., 7b in Eqs. (|T4|) and 
(fT5|) now has the same dependence on Rabi frequency as A. This is why ares for the 
resonances in Fig. [5] is independent of Rabi frequency. However, as shown in Fig. [3l a 
significant change in the bound state wavefunction can lead to a strong change in a^es- 

4-2. ^'^Rb: bound-bound coupling 

^''Rb has been used in a wide range of studies in ultracold gases. A large number 
of resonances have been observed in this system [53], with many studies utilising the 
comparatively wide (A = 0.21 G) resonance in the a + a channel at Bq = 1007 G |54j . 
There are a number of resonances grouped close to each other around 9G and 18 G, 
in channels corresponding to the zero field (/i = 1) + (/2 = 2) limit P, [55l [56l [32] . 
These resonances were used in the demonstration and theoretical analysis of rf-dressed 
Feshbach resonances in Ref. [1]. We also note the calculations presented in Ref. [23]. In 
Ref. [I], the primary effect of the rf was to couple the bound states causing each of the 
nearby resonances. Bound-bound coupling has a substantially larger Franck-Condon 
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Figure 5. Real part of the scattering length (upper panel) and two-body loss rate 
coefScient (lower panel) as a function of magnetic field. Here, we use the ja + fo, iVo) 
entrance channel and polarised rf with an oscillation frequency of 77.5 MHz to create 
a resonance near 570 G. The Rabi frequencies used are given in MHz. The width of 
the loss feature increases strongly with Rabi frequency, as does the background loss 
created by rf coupling to lower exit channels. However, the maximum available change 
in the scattering length remains constant. 



overlap than bound-free coupling, and so significant control of scattering properties 
can be achieved with correspondingly less rf power. For the calculations shown in 
this section, the Rabi frequency on the a ^ h transition is 38 kHz, corresponding to 



All collision channels of ^''Rb have a similar value of Obg, because of the similarity 
of the singlet and triplet scattering lengths. Another consequence of this is that the 
highest vibrational bound states for all possible F have approximately the same binding 
energy at i? = 0. Figure |6^ shows the F = 1,2, and 3 bound states closest to the 
(/i = 1) + (/2 = 2) threshold. At zero field, the F = 1 and F = 3 bound states are 
degenerate, with the F = 2 state approximately h x 1 MHz deeper. The atomic and 
molecular states are split into their Zeeman components at nonzero field. At magnetic 
fields above 1 G, each bound state is strongly concentrated (> 90 %) in a single channel 
a + f3, and is labelled (a/3). The energy of each bound state then remains at a constant 
offset from the threshold energy of its corresponding channel. Also, the Zeeman effect 
is close to linear within the range of magnetic field considered here. Consequently, for a 
fixed magnetic field, each pair of adjacent channels or bound states (of the same /i and 
/2) have approximately the same energy gap. This is illustrated in Fig. [6]d, which shows 
the bound states of the a + d, a + e and a + / channels. Here, energies are shown relative 
to the threshold of each channel. The three spectra are nearly identical, highlighting 
the symmetries inherent in the Zeeman interaction. The combinations of atomic states 
from which channels of each Mt can be formed, and the bound state energies at i? = 0, 
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Figure 6. Bound state structure and Feshbach resonances of ® Rb. Energies are 
shown as a function of magnetic field. Bound states are shown as sohd lines, and 
labelled with the channel (a/3) in which they are concentrated. Collision thresholds 
are shown with dotted lines. In (a), energies are shown relative to the a + e entrance 
channel threshold. Resonances occur when a bound state crosses a collision threshold 
of the same Mt- These are shown with points and labelled with the entrance channel 
a + (3. Note that we only show the thresholds of channels containing a resonance. 
Colors indicate different Mt, as labelled. The Zeeman effect is close to linear in the 
field range shown, and each bound state is concentrated in a single channel. This gives 
rise to the grouping of resonances around 9 G and 18 G. In (b), the bound state spectra 
of the a + d, a + e and a + f channels are shown relative to their own thresholds. In 
(c), the dressed channel energies are shown relative to that of the |a + e, Nq) threshold. 
Here, the rf frequency is 5.5 MHz. The bound state structure shown in (a) and (b) 
results in a number of bound states becoming degenerate with each other and the 
\a + e, A^o) threshold near 9 G. The bound state causing the a + e resonance is coupled 
most strongly to states in the \a + f,NozL 1) channels, which are indicated with the 
dashed circle in panel (c). 

explain the existence of the three Feshbach resonances near 9 G and the eight near 18 G 
shown in Fig. [6^ [H |57] . 

Another consequence of the Zeeman effect being predominantly linear is that several 
channels can be resonantly coupled by a single rf frequency, as shown by the avoided 
crossings in Fig. O It is also possible to couple several bound states, and use the 
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magnetic field to tune these coupled states through a collision threshold. We consider 
atoms colliding in the |a + e, A^o) entrance channel, and choose an rf frequency close to 
the atomic Zeeman splitting at 9 G. The resulting energies of the \a + f,No + 1) and 
|a + /, A'o — 1) bound states relative to the \a + e, Nq) threshold are shown in Fig. 
for an rf frequency of a; = 27r x 5.5 MHz. There is near-degeneracy between the bound 
states causing the 9 G resonances in the \a + e, Nq) and \a + f,NQ — 1) channels, and 
that causing the 18 G \a + f,NQ + 1) resonance. Consequently, a pair colliding in the 
entrance channel are linked by spin exchange to a set of strongly rf-coupled bound states. 

Several other bound states contribute to the rf-dressed scattering around 9 G, 
leading to complicated variation of scattering properties with magnetic field and rf 
frequency. These are not shown in Fig. ^ for reasons of clarity. The real part of the 
scattering length a{B,u), and the two-body loss coefficient K2{B,u), are plotted in 
Fig. [71 The broad, horizontal band in both plots corresponds to the undressed a + e 
resonance. Coupling of the underlying bound state to others creates several secondary 
features at magnetic fields that depend on the rf frequency. When this field is close 
to 9.1 G, the a + e resonance is split. The rf-induced avoided crossing then changes 
the magnetic field at which each bound state crosses the entrance channel threshold. 
Although it is therefore possible to suppress losses at the centre of a resonance or move 
its location, this will in general also suppress or move the resonant enhancement of 
the scattering length. This is shown by the variation in scattering length of Fig. [7^, 
which follows the same trends as the losses in Fig. [7)d. Loss rates generally have a peak 
close to Feshbach resonances, and are therefore commonly used for deducing resonance 
locations. In Fig. [7)d, we reproduce the experimental data of Ref. [1], showing its close 
agreement with our calculations. Each point represents the centre of a loss feature. The 
bound state giving rise to each feature is identified by the rf-dressed channel |(aia2), N) 
in which it is concentrated. The undressed, magnetically tunable resonances produced 
by each of these bound states are shown in Fig. [6^. 

For the Rabi frequency of 38 kHz used here, the bound-free coupling induced by the 
rf is minimal, and bound-bound coupling dominates due to the larger Franck-Condon 
overlap. Because of the dominance of bound-bound coupling, the locations of the rf 
dressed resonances can be reproduced well using a simplified Hamiltonian consisting 
only of the bound states which are resonantly coupled near 9 G. For this, we assume 
that the bound state energies are given simply by the undressed energy shifted by 
the number of photons that makes it near-resonant with the \a + e,No) threshold, e.g. 
E\{bg),N) = -E'(f)g) + (A^ — Nq)Jtjjj. Thls neglects any shifts in the threshold energies due 
to the intensity of the light. We also use the atomic Rabi frequencies to represent the 
coupling between bound states, e.g. fi(bg)^(cg) = ^^b^c- This neglects any differences 
in the molecular wavefunctions, the overlap of which should be close to unity for this 
case of narrow resonances. We diagonalise the Hamiltonian and find the fields and rf 
frequencies for which the dressed bound states become degenerate with the \a + e,N) 
threshold. These are shown as dashed lines in Fig. [71 The good agreement of this 
simpler calculation confirms the dominance of bound-bound coupling in the effects of 
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Figure 7. (a) Real part of the scattering length and (b) two-body loss coefBcicnt. K2, 
as a function of magnetic field and rf frequency. Black circles show the experimental 
data of Ref . [I] , which identified the rf-dressed resonances by the centre of loss features 
corresponding to the maximum of 7^2- Dashed lines correspond to the locations 
predicted using the simple model of Sec. 14.21 which includes only bound states. The 
band of sudden changes in the middle of each figure corresponds to the region in which 
the resonances of atomic transitions are crossed. This makes both calculations and 
measurements impractical. 



the rf on the scattering properties. 

We note that unhke the losses found for ^Li in Sec. 14.11 which are into exit 
channels coupled by the rf, the losses here are due to the properties of the magnetically 
tunable resonances themselves. All of these resonances are closed-channel dominated 
and strongly decaying, due to inelastic spin relaxation loss into lower channels of the 
same Mt [1]. Each has a width of order niG and a resonance length of ares ^ abg/2. 
Consequently, the strongest losses in Fig. [7] are independent of Rabi frequency, and a 
secondary feature due to a bound state from a channel with photons grows in strength 
proportional to fi^l^"^"'. 

4.3. ^Li: free-free coupling 

It is possible to control scattering properties by creating rf-dressed atomic states. This is 
done routinely in the creation of clock states [12]. For rf of a frequency close to an atomic 
transition, the dressed atoms have a substantial admixture in the two coupled Zeeman 
states. A collision between two dressed atoms then samples the scattering properties of 
all the corresponding undressed two-body channels. For such free-free coupling to have 
an effect requires the scattering length in each undressed channel to be substantially 
different. In ^^Rb, for example, all channels have very similar background scattering 
lengths. This makes the collision of two atoms in dressed states similar to that of two 
atoms in the absence of rf. By contrast, in Fig. |8]we show the example of ^Li scattering. 
The broad, overlapping Feshbach resonances in the a + b, a + c and b + c channels make 
the difference in scattering lengths between the channels have a nontrivial variation with 
magnetic field. Large changes in scattering length are produced for rf frequencies close 
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Figure 8. Scattering length and loss rate as a function of rf frequency uj/2tt. The 
magnetic field is 250 G, and the cr~ rf has a Rabi frequency of 1 MHz. Significant 
changes in scattering properties occur close to the resonant frequencies of the a -H- 6 
and 6 o c atomic transitions. A two-photon a -f-)- c transition is also driven. The 
loss properties of each channel change dramatically as the transition frequencies are 
crossed. Away from the transition frequencies, small changes in scattering length are 
achievable with low loss rates. 

to the atomic transitions, over a range comparable to the Rabi frequency. 

The scheme described in this subsection does not create a Feshbach resonance. It 
does, however, sample the scattering lengths in different Zeeman channels, including 
any Feshbach resonances those channels support. Very close to the atomic transition 
frequencies, losses will be a problem for this scheme, as shown in the lower panel of Fig. |8l 
However, such a method could be useful for creating small changes in the scattering 
length, and may remain an option in the absence of bound states suitable for using the 
methods of the previous two subsections. In the example given here, a 5% change in 
scattering length can be achieved with loss rates of order 10^^'^ cm^/s. 

5. Conclusions 

We have developed a technique for studying the collisions of cold atoms in the presence 
of radio-frequency radiation. Building on the three-parameter model of Feshbach 
resonances presented in Ref. [33], we have incorporated a frame transformation to a 
basis in which Zeeman states are coupled together by rf radiation. Our studies have 
shown that rf provides useful control capabilities when used in conjunction with one or 
more Feshbach resonances - either by coupling together bound states with which the 
colliding atoms interact [1], or by directly coupling the colliding pair to a molecular 
state. For the latter case, we have found that the use of a halo molecule reduces the rf 
power required for control. Also, some ranges of rf frequency and magnetic field exist 
for which a bound state can be coupled without causing losses. This requires the use of 
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circularly polarised rf. 

The accuracy and speed with which rf can be controlled, and the potential of tuning 
several scattering lengths of a multicomponent gas, give our work relevance to many 
current experimental programs in atomic and molecular collisions. Other candidates for 
such an enhancement of control possibilities include ^"K, which has broad resonances in 
the a + b and a + c channels. Another example more like the ®^Rb calculations of Sec. 14.21 
is ^Li-^^K, which has several narrow resonances around 160 G. The a + a channel could 
allow bound-bound coupling without the strong decay inherent to the excited state *^Rb 
resonances. 
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